Electron-Electron Interactions on the Edge States of Graphene: A Many Body 

Configuration Interaction Study 
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We have studied zigzag and armchair graphene nano ribbons (GNRs), described by the Hubbard 
Hamiltonian using quantum many body configuration interaction methods. Due to finite termina- 
tion, we find that the bipartite nature of the graphene lattice gets destroyed at the edges making 
the ground state of the zigzag GNRs a high spin state, whereas the ground state of the armchair 
GNRs remains a singlet. Our calculations of charge and spin densities suggest that, although the 
electron density prefers to accumulate on the edges, instead of spin polarization, the up and down 
spins prefer to mix throughout the GNR lattice. While the many body charge gap results in in- 
sulating behavior for both kinds of GNRs, the conduction upon application of electric field is still 
possible through the edge channels because of their high electron density. Analysis of optical states 
suggest differences in quantum efficiency of luminescence for zigzag and armchair GNRs, which can 
be probed by simple experiments. 
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Nanomatcrials of carbon of different dimensionalities 
have been a subject of interest over the past few decades 
due to their potential applications in various nanoscale 
electronic devices 1 ^! 3 -. Compared to fullerenes and car- 
bon nanotubes, which are effectively zero and one di- 
mensional in nature, the two dimensional flat mono- 
layer of carbon atoms packed into a honeycomb lat- 
tice, i.e., graphene, has started gaining prominence very 
recently owing to the recent progress in experimen- 
tal techniques 4 ^^. Because of its sophisticated low- 
dimensional electronic properties and huge application 
possibility, it has attracted a big scientific army to ex- 
plore it in various aspects ^ 10 ' 11 . 

The size and geometry of the nanoscale carbon sys- 
tems govern their electronic properties. Recent progress 
in experiments allows to make finite size graphene 
layer, termed as graphene nanoribbons (GNRs) with 
varying widths, either by cutting mechanically exfoli- 
ated graphenes 4 ^!^ and patterning by electron beam 
lithography^, or by controlling the epitaxial growth of 
graphene a 13 i 14 . Different possibilities of geometrical ter- 
mination of graphene layer give rise to two different edge 
geometries, namely, zigzag and armchair edges, differ- 
ing largely in their electronic properties. These different 
edge geometries have been modeled by imposing different 
boundary conditions on Schrodinger's equation within 
the tight-binding limi t 15 ' 16 i 17 i 18 or on the Dirac equation 
for two-dimensional massless fermions with an effective 
speed of ligh t 19 ' 20 i 21 ' 22 in previous studies. These also 
have been extensively studied using density functional 
theor y 23 i 24 . However, the proper many-body description 
of the GNR is still lacking. To obtain new inroads in the 
rich physics of one atom thick finite graphene material, 
for the first time, we use configuration-interaction (CI) 
method based on many-electron theory with proper in- 
clusion of electron interactions term which plays crucial 
role in low-dimensional systems. 

Our CI approach with correlation parameters in 




FIG. 1: The unit cell of GNR. The translations along x and 
y-axis give the zigzag and armchair edge GNRs. The num- 
ber of atoms at the zigzag terminal of an armchair GNR are 
represented by the integers. 



Hamiltonian obtained semiempirically, can directly con- 
nect to the real description of the system under inves- 
tigation. We study both the zigzag and armchair edges 
with a fixed width and vary their size by translating it 
in one direction up to a large limit to reduce finite size 
effects. Since, the full many-body CI calculation is not 
possible for very large system, we map the whole system 
into a complete active space (CAS) which captures the 
low energy states accurately. We have also considered 
the Single CI (SCI) which shows the same qualitative 
features as obtained from CAS-CI. 

We employ Hubbard Hamiltonian for the GNRs, 



+ U ^n^mi (1) 



with standard notations. The nearest neighbor hopping 
integral ty is considered to be 2.4 eV which is standard 
for the C-C bond. We consider U as the difference be- 
tween the first ionization energy and the electron affinity 
for carbon, which comes out to be ~ 9.66 eV. Due to 
finite termination of the GNRs. the chemical nature of 
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FIG. 2: The lowest gap in tight binding calculations as a 
function of inverse system size {1/N) for both zigzag (circle) 
and armchair (square) GNRs with e e d ge equals to (a) — 2.0 
and (6)0. 



TABLE I: The ground state energies (E(S Z )) per atom for 
two different system sizes (N) for both zigzag and armchair 
GNRs at half filling with different e e d ge (eV) values. Since all 
total spin states (5*°') have projection into 5* ot < S tot , the 
ground state spin is determined by the highest Sl ot value with 
ground state energy. 
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the edge atoms changes abruptly because of unsatisfied 
valence of carbon atoms. As a result, the lattice deviates 
from its bipartite nature at the edges. So, to differentiate 
these edge atoms from the bulk, we model our system by 
putting non zero negative onsite energy, e e d ge at the edge 
atoms to mimic hydrogen passivation and ej = for the 
bulk. Although, the onsite energies at edges are expected 
to be modified compared to the bulk as a difference be- 
tween the electronegativities of carbon and hydrogen^, 
they have been neglected in most of the studies so far. 
Moreover, the scanning tunneling microscopy (STM) im- 
ages of graphite shows bright stripes along the edges, 
suggesting more electron density on the edge atoms com- 
pared to the bulk2£. This observation motivates us to 
simulate the GNRs with negative onsite energies at the 
edge carbon atoms to create a potential well which can 
trap the electrons on the edges. We compare and con- 
trast our model with the bipartite lattice model, where 
all the lattice points are assigned = 0, both at the tight 
binding (TB) and at the CI level of calculations. 

Fig.l displays the GNRs under investigation. The 
translations along x-axis and y-axis produce zigzag and 
armchair GNRs of different sizes respectively. We con- 
sider the armchair GNRs with 3p number of atoms 
constituting the two zigzag terminals, where p (= 3 
in our case) is an integer. In Fig. 2, we plot the TB 
(U = 0) gap for both types of GNRs with and with- 
out consideration of edge hydrogen passivation energies, 
with inverse system size. The calculated TB gap re- 
produces the well established previous observations of 
zigzag being metallic and armchair (3p) being semicon- 
ducting GNRs^^^^^^*^. Howeverj inclu _ 

sion of electron correlations are expected to open up a 
gap reducing the kinetic stabilization. 

Based on the TB wave functions, we perform the CI 
calculations over the electron correlations within many 
body formalism. The z-component of total spin (5'°') is 
used as a quantum number together with the N number 
of electrons. We vary S 1 * * from zero to higher values for 
half-filled systems. In CAS-CI, we consider all possible 
configurations within a small energy window (the active 
space) and obtain the CI matrix of order 4900. We have 
verified our results with 63504 states in certain cases. 



We have also considered singles CI (SCI) calculations by 
varying the number of occupied and unoccupied single 
particle levels. However, since the CAS-CI includes all 
the configurations within a certain energy cut-off, the re- 
sulting low-energy states are more size consistent than 
those from the SCI and unless otherwise stated, all the 
results reported below are obtained using CAS-CI ap- 
proach. For any kind of energy cut off with S 1 "* — 0, we 
find that the diagonals of the CI matrix appears in the 
order of after correction for the levels beyond the cut 
off. 

Our many-body calculations for half-filled GNRs with 
different choices of e e d ge value and with varying S z states, 
lead to results, which have been anticipate d 23 ! 27 but not 
observed till now: that, the ground state of the zigzag 
GNR is always a high spin state, whereas, it remains in 
a singlet state for the armchair GNRs. In Tablel, the 
numbers in parenthesis for various size GNRs are the 
total ground state spins. This observation remains con- 
sistent for all nonzero negative e e d ge values. However, 
for e e d ge = 0, both the GNRs show singlet ground state, 
which is expected when the lattice is bipartite with same 
number of atoms in the two sublattices 2 ^. These obser- 
vations are consistent with the gradual increase of the 
system size for a fixed width. To correlate the exper- 
imental observation^, we calculate the charge density 
over all the atoms for both zigzag and armchair GNRs 
with N = 104, with e edge = and -2.0 (Fig.3). The 
tedge = value shows almost same electron densities on 
all the atoms, whereas non-zero e e dge leads to fairly large 
charge densities on the edge atoms (maxima in Fig.3) for 
both types of GNRs as can be seen from Fig.3. This is 
the result of unsatisfied coordination of the edge atoms, 
passivated by hydrogen. Interestingly, the charge accu- 
mulation at the zigzag edges are more than those at the 
armchair edges. Note that, the minima in the charge 
density plots correspond to the atoms directly connected 
to the edge atoms. This reduction in charge densities on 
the left and right neighbors of edge atoms gives kinetic 
stability to the extra electrons accumulated at the edge 
atoms. 

To understand the nature of the spin at the edges, 
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FIG. 3: Charge density on each and every site for both zigzag 
(a) and armchair (b) GNRs with e e d ge = —2.0 (grey lines with 
circles) and (solid lines without any symbol) on the edge 
atoms of the system with iV = 104. Insets present the contour 
plot of the charge density for the corresponding systems with 

tedge — 2.0. 

we calculate the spin density over all the atoms, which 
suggests that both zigzag and armchair GNRs do not 
prefer any dominant spin at the edges as well as on the 
bulk atoms for both zero and non-zero e e dge values. In- 
stead, a proper many body consideration shows that, the 
GNRs prefer to have a mixture of both up and down 
spin density at the edges. This observation clearly con- 
tradicts previous density functional studies which pre- 
dict the dominance of up spins and down spins on cither 
edges of zigzag GNR a 23 i 24 . Those studies conclude this 
observation as the property of bipartite lattice with two 
different sublattice points. However, the presence of pas- 
sivating hydrogen atoms on the edges would destroy the 
bipartite nature of the lattice which, however, have been 
completely ignored. Moreover, the one electron defini- 
tion of the exchange correlations in DFT studies cannot 
capture the dominant electron electron interactions in 
GNRs. The more realistic many body calculations infer 
that, although the electrons have the tendency to accu- 
mulate at the edges, a net spin polarization of the edges 
is highly improbable. 

Now to investigate the possibility of electron conduc- 
tion in the different GNRs, we calculate the charge gap. 
The many-body charge excitation gap is defined as the 
difference between the energy required to add (/i+) and 
remove (/i-) electrons from the ground stated, 

^■charge = M+ — (2) 

where ^+ = E(N+1)-E(N) and /x_ = E(N)-E(N-l). 
E(N), E(N + 1) and E(N - 1) arc the energies of the 
half-filled system and the systems with one extra and one 



less electron respectively. The charge gap results from 
our calculations show Mott insulating behavior for both 
zigzag and armchair GNRs even for very large system 
sizes with more than thousand atoms as observed in ear- 
lier many body study^i. This result is expected since the 
kinetic stabilization is suppressed by strong electron cor- 
relations, which is unexpected within one-electron theo- 
ries. For a half-filled system, with large electron densities 
at the edges, the conduction, however, is still possible 
through edge channels. 

Motivated by the difference in the spin multiplicities 
of the ground and low-lying excited states between the 
zigzag and armchair GNRs, we study their optical prop- 
erties. We have calculated the transverse and longitudi- 
nal transition dipole moments, /zt and /il respectively, 
for the excitations from the ground state to the opti- 
cally allowed excited states and present /i 2 = + n\ 
for N = 104 in Table 2. As can be seen from Table 
1, the energy levels in zigzag GNRs are more closely 
spaced than those in the armchair GNRs. Similar to 
the ground state, the first few excited states of zigzag 
GNRs are also of higher spin states. However, due to 
nonzero e ee i ge , the high spin ground states are not de- 
generate with different 5 , ' ot components. With nonzero 
gap values, these low energy excitations show large ab- 
sorption cross sections. Moreover, the optically allowed 
emissions from the excited states are expected to show 
very high quantum efficiency (QE) in the luminescence 
spectroscopy, since according to Kasha's rule 3 ^ the spin 
allowed optical emission can occur from high spin excited 
state to the high spin ground state. In the case of arm- 
chair GNRs, however, the singlet ground state is followed 
by a few high spin excited states. This dipole forbidden 
excited states below the allowed optical excitation orig- 
inate from the strong electron electron interaction s 31 ' 32 . 
Thus, it is expected that, the optically excited state in 
armchair GNRs would decay to the low energy dipole 
forbidden states through magnon emission, thereby pre- 
venting the radiative transition considerably. As a result, 
the armchair GNRs are expected to show very low QE in 
luminescence spectroscopy. These results are all consis- 
tent with the gradual increase of system size with fixed 
width. Thus, from our study on optical properties, we 
are able to propose a simple experimental technique of 
measuring luminescence to differentiate the zigzag and 
armchair GNRs. 

While our discussions so far are based on short range 
Hubbard Hamiltonian, it is known that the one dimen- 
sional conjugated carbon systems are best described by 
long range Coulomb interactions 3 ^. For graphenc class of 
systems, however, even the on-site Hubbard correlations 
have hardly been considerecU^^^^^^ 2 ^ 3 * 2 ^. 
There are also suggestions that the Hubbard repulsion 
can at best be 3 — AeV in graphene^ 3 ., three times smaller 
in magnitude than what is considered for the conven- 
tional conjugated polyemcrs^ 2 -. Interestingly, we have 
verified that all our observations described earlier are 
consistent even with U = 3eV. However, to test whether 
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TABLE II: The square of the transition dipole moment (/i 2 ) 
for transitions from the ground state to the few optically al- 
lowed excited states with the corresponding excitation gaps 
for the system with N = 104 in both zigzag and armchair 
GNRs with e edge = -2.0 
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Gap (el/) 




0.016 


2.801 




0.906 


3.770 


0.071 


0.440 




0.969 


0.001 


0.073 


3.069 




1.060 


3.027 


0.115 


0.168 




1.136 


0.002 



longer range interactions exist in such systems, we have 
investigated these GNRs with long-range Coulombic in- 
teractions using PPP Hamiltonian within two decay pro- 
file schemes originally given by Ohno and by Mataga- 
Nishimoto^i. Our results indicate that for both these 
potential profiles, the ground state has large contribution 
from either high-energy singles configuration (in case of 
SCI) or configurations with two or three electron exci- 
tations within the energy window for the CAS-CI calcu- 
lations. Contributions from these configurations in fact 
raise the energy of the state with finite magnetic mo- 
ment. This hence results in both zigzag and armchair 



GNRs having a singlet ground state, while the lowest 
excitations become magnetic. 

In summary, the quantum many-body configuration 
interaction method captures the low-energy properties of 
the nano scale systems like graphene and our develop- 
ment of CI method allows us to handle fairly large sys- 
tems which can effectively be considered as infinite lat- 
tice. Within the Hubbard model, the ground state of the 
zigzag GNRs is a nonzero spin state, while for armchair 
GNRs it remains a singlet. Though the hydrogen passi- 
vation of the edge atoms leads to higher charge density 
at the edges, the electronic correlations mix the up spin 
and down spin throughout the GNR lattices, instead of 
making the edges spin polarized. We propose that, the 
zigzag GNRs can be differentiated from armchair GNRs 
from the higher quantum efficiency of luminescence orig- 
inating from their magnetic ground and low lying excited 
states. It would also give an indication of the importance 
of the long range Coulomb interactions in these class of 
systems. Our findings on GNRs suggest rich low-energy 
physics and provoke further studies within many body 
limit for these class of low-dimensional systems. 
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